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f—^ , In this paper a recently developed projector-based renormalization method (PRM) for many- 

^— V ■ particle Hamiltonians is applied to the periodic Anderson model (PAM) with the aim to describe 

fvj ' heavy Fermion behavior. In this method high-energetic excitation operators instead of high energetic 

states are eliminated. We arrive at an effective Hamiltonian for a quasi-free system which consists 
1-" . of two non-interacting heavy-quasiparticle bands. The resulting renormalization equations for the 

~T ' parameters of the Hamiltonian are valid for large as well as small degeneracy Vf of the angular 

momentum. An expansion in l/vf is avoided. Within an additional approximation which adapts 
f~^ , the idea of a fixed renormalized / level e/, we obtain coupled equations for if and the averaged 

/ occupation (ny). These equations resemble to a certain extent those of the usual slave boson 

mean-field (SB) treatment. In particular, for large Vf the results for the PRM and the SB approach 

agree perfectly whereas considerable differences are found for small Vf. 
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Y ■ I. INTRODUCTION 

^'. 

C ' In comparison to ordinary metals metallic heavy fermion systems have remarkable low-temperature properties [ij: 

both the conduction electron specific heat and the magnetic susceptibility can be two or more orders of magnitude 
larger than in normal metals though the ratio of both quantities is similar to that of usual metals. Usually, in 
metals an increasing resistivity p{T) with increasing temperature is observed. In contrast, a much richer behavior 
\^ ' is found in heavy fermion systems: At higher temperatures p{T) only changes slightly and might even increase 
with decreasing temperature. Below a characteristic coherence temperature a strong decrease of the resistivity with 
decreasing temperatures is observed. At very low temperatures, a T^ dependence of the temperature is found. Another 
\r^ ' important finding is that a correspondence between the low-energy excitations of heavy fermion systems and those of 
\^ , a free electron gas with properly renormalized parameters can be established. Therefore, the high density of states at 
^^ the Fermi surface observed in heavy fermion systems implies an effective mass of the (heavy) quasiparticles which is 
\l , some hundred times larger than the free electron mass. 

^^ ■ Prototype heavy fermion systems like CeAls and UPts contain rare-earth or actinide elements. Thus, the basic 

microscopic model for the investigation of such materials is believed to be the periodic Anderson model (PAM) which 
describes the interaction between nearly localized, strongly correlated/ electrons and conduction electrons Q- Within 
a simplified version the Hamiltonian of the PAM can be written as 

"§ ■ n = Ho+Hi, (1) 
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>< : Hi = ^ ^ Fk (/Lck™ e^'-^' + h.c. 



Here, i is the 4/ or 5/ site index, k is the conduction electron wave vector, and Vk is the hybridization matrix 
element between conduction and localized electrons, e/ and Sk, both measured from the chemical potential /i, are 
the excitation energies for / and conduction electrons, respectively. As a simplification, both types of electrons are 
assumed to have the same angular momentum index m with i/f values, m — \...Vj. Finally, the local Coulomb 
repulsion Uf at / sites has been assumed to be infinitely large so that localized sites can either be empty or singly 

occupied, i.e., the Hubbard operators /j„j are defined by 

/L=/L n (i-"L) (2) 



where nf^ = fj^fim- The unexpected and exciting properties of the PAM ^ are mainly due to the presence of 
the strong correlations at / sites. In turn the strong correlations also cause the great difficulties in any theoretical 
treatment of the model. In the present approach the correlations are taken care of by the Hubbard operators (0 
which do not obey the usual fermionic anticommutator relations. Instead one has 

where 
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The quantity Vim can be interpreted as a local projection operator at / site i on f states which are either empty 
or singly occupied with one electron with index m. Also it is helpful to introduce separately the projection operator 
7^o(*) on the empty / state at site i and the projection operator fij^ on the singly occupied / state when one electron 
with index rn is present. Vim can be rewritten as 

V„n = Po(*)+nL = 1- E ^irn (4) 



where we have defined 



'P,{^) = Y[(^-.fLf^r~n), (5) 
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The second equation in Q is the completeness relation for / electrons at site i. 

For the case of vanishing Coulomb repulsion Uf the PAM ^ is equivalent to the Fano- Anderson model y, |j| 
which can be easily solved (see, for example, appendix A). However, much of the physics of the correlated model 
can also be understood in terms of a renormalization of the parameters of the uncorrelated Fano- Anderson model. 
Various theoretical methods have been developed in the past to generate such renormalized Hamiltonians, for instance 
the Gutzwiller projection [y or the slave-boson mean-field (SB) theory [Q, Q- Here we use a recently developed 
[3 projector-based renormalization method (PRM) to map the PAM to a free system consisting of two bands of 
uncorrelated quasi-particles. Furthermore, we avoid an expansion with respect to the degeneracy Vf of the angular 
momentum and take all l/i^/ corrections into account. 

The PRM has already been applied before to the PAM in Ref. &. However, in the present approach the treatment 
from Ref. |3 will be improved in various points: (i) The PRM is performed in a completely non-perturbative manner, 
(ii) All 1/vf corrections are taken into account, (iii) The dispersion of both quasiparticle bands is considered. 

Furthermore, we shall compare the results of the PRM with those of the SB treatment in much more detail. 

The paper is organized as follows. First, in Sec. |n]we briefly repeat the recently developed PRM 8]. In Sec. IIIII 
the PRM is applied to the PAM whereby the renormalization equations for the model parameters are derived non- 
perturbatively. An analytical solution of the renormalization equations is found in Sec. IIVI using a constant renormal- 
ized / level if. Furthermore, we compare our results with the solutions of the SB theory. Finally, our conclusions are 
presented in Sec. |V| 

II. PROJECTOR-BASED RENORMALIZATION METHOD (PRM) 

The PRM Q starts from a decomposition of a given many-particle Hamiltonian Ti into an unperturbed part Tig 
and into a perturbation Tii, 

n = Ho+Hi. (7) 

We assume that the eigenvalue problem of TYq is known 

no\n) = i?(°)|n). (8) 



Til is the interaction. Its presence usually prevents the exact solution of the eigenvalue problem of the full Hamiltonian. 
Let us define a projection operator P> by 

Pa-4 = Y^ \n){m\{n\A\m). (9) 
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Note that P^ is a super-operator acting on usual operators A of the unitary space. It projects on those parts of A 
which are formed by all dyads |'^)(ti| with energy differences \En — Em \ less or equal to a given cutoff A, where A is 
smaller than the cutoff A of the original model. Note that in neither \n) nor \m) have to be low-energy eigenstates 
of Ho- However, their energy difference has to be restricted to values < A. Furthermore, it is useful to define the 
projection operator 

Qa = 1-Pa (10) 

which is orthogonal to Pa. Qa projects on high energy transitions larger than the cutoff A. 

The goal of the present method is to transform the initial Hamiltonian Ti. (with a large energy cutoff A) into an 
effective Hamiltonian Tix which has no matrix elements belonging to transitions larger than A. This is achieved by an 
unitary transformation so that the effective Hamiltonian will have the same eigenspectrum as the original Hamiltonian 
Ti. However, as it will turn out, the method is especially suitable to describe low-energy excitations of the system. 
H\ is defined by 

Hx = e^^ne-^\ (11) 

The generator X\ of the transformation has to be anti-Hermitian, Xj^ — ~X\, so that Tix is Hermitian for any A. 
We look for an appropriate generator A"a so that Tix has no matrix elements belonging to transitions larger than A. 
This means that the following condition 

QaHa = (12) 

has to be fulfilled. Eq. (|12|) will be used below to specify Xx- In contrast to Ref. [3], where Tix was evaluated 
perturbatively, the transformation (|ll|l will be treated non-perturbatively. 

Next we discuss the elimination procedure for the interaction Tii. Instead of transforming the Hamiltonian in 
one step as in Eq. Hll() the transformation will be done successively. Or more formally spoken, instead of applying 
the elimination of high-energy excitations in one step a sequence of stepwise transformations is used in order to 
obtain an effectively diagonal model. This procedure resembles Wegner's fiow equation method ^ and the similarity 
renormalization 10] in some aspects. In the PRM approach difference equations for the A dependence of the parameters 
of the Hamiltonian are derived. They will be called renormalization equations. To find these equations we start from 
the renormalized Hamiltonian 

"Hx = Wo,A+Wl,A (13) 

after all excitations with energy differences larger than A have been eliminated. Next we consider an additional 
renormalization of TLx by eliminating all excitations inside an energy shell between A and a smaller energy cutoff 
(A — AA) where AA > 0. The new Hamiltonian Ti.(x-AX) is obtained by an unitary transformation similar to that of 
Eq. mi) 

H(A-AA) = e^^-^^ T^A e-^^'^^ (14) 

where ATa.aa is determined by 

Q(A-AA)'^(A-AA) = 0. (15) 

Note that there are two strategies to exploit Eq. H15|) in order to determine the generator Xa,aa of the unitary 
transformation H14f) . The most straightforward route is to analyze Eqs. (|14|l and (|15|) in perturbation theory as it was 
done in Refs.lSl and I 111 Here, we want to perform the renormalization step from A to (A — AA) in a non-perturbative 
way. 

Eqs. H14|l and H15|) describe the renormalization of the Hamiltonian by decreasing the cutoff from A to (A — AA) and 
can be used to derive difference equations for the A-dependence of the Hamiltonian. The resulting equations for the 
parameters of the Hamiltonian will be called renormalization equations. Their solution depends on the initial values 
of the parameters of the Hamiltonian and fixes the final Hamiltonian in the limit A ^- 0. Note that for A ^ the 
resulting Hamiltonian only consists of the unperturbed part TIq rx^o}- The interaction Tiiix^o) vanishes since it is 
completely used up in the renormalization procedure. Thus, an effectively diagonal Hamiltonian is obtained. 



III. RENORMALIZATION OF THE PERIODIC ANDERSON MODEL (PAM) 

The PRM described above will be applied in this section to the PAM iQJ. As an illustration of the method the 
Fano- Anderson model is discussed as a further application of the PRM in appendix EI This model can be considered 
as a PAM without electronic correlations. It turns out that the renormalization of the full PAM |^ is somewhat 
similar to that of the uncorrelated model. However, in the Anderson- Fano model the elimination of excitations with 
energies larger than A can be done in one step. For the PAM (^ the / electron one-particle energy Sf will also 
be renormalized. This is due to the presence of strong local correlations at / sites in the PAM Q). Therefore, the 
elimination procedure has to be done stepwise by repeatedly integrating over small energy steps of width AA. In this 
way one is led to renormalization equations for the parameters of the model in terms of difference equations which 
have to be solved. 



A. Renormalization ansatz 

Let us start by formally writing down the effective Hamiltonian Hx = e^^T-ie~^^ = 7io,A + 'Wi,a for the periodic 
Anderson model after all excitations with energy differences larger than A have been eliminated. By comparing with 
the starting model ^ one might be in favor of choosing an unperturbed part 7io,A which contains correlated /- 
electrons /J^, fim as in Eq. (^. However, because the eigenvalue problem of such an Hamiltonian would not exactly 
be solvable, we prefer to start from an uncorrelated Hamiltonian where only usual Fermi operators fj^ , fim enter but 
keep the correlations in the renormalized interaction Tii^x, i.e. 



T^x = Hqx +Hi, 



where 



'Ho,x = e/,A XI /km/km + X! V^ (/L/km) ^^^ + XI '^k,A cLckm + £^A, (16) 

k,/n k,m k.TTi 

Wi,A - PaWi - X^kPA (/Lckm +h.c.) . (17) 



Here /ij,„, /j[„ are Fourier transformed / operators. 
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Moreover, in Hlt)|) we have discriminated between local, 

[jmJmJj^ '■= 'TJ / , /km -/km = T7 /_^ Jim i*™ ' (19) 



k 



and nonlocal, 



l^/km/kmj^^ :— ^ 2^ JimJjm^ ' ^ — JkrnJ^"^ [jmJmJi^, (20) 

/ particle-hole excitations in order to properly take into account the strong Coulomb interaction at local / sites. 
Due to renormalization processes the one-particle energies Sf^x and ek,A in Eq- H16() depend on the cutoff energy 
A. Moreover, two new parameters enter: Ak,A describes the / dispersion due to the hybridization of / electrons at 
different sites i ^ j, and Ex is an additional energy shift. Finally, the projector Pa in (|17|l guarantees that only 
excitations contribute to Tii^x which have energies (with respect to Ti-o^x) which are smaller than A. The initial 
parameter values of the original model (at A = A) are 

£/,(A=A) = £/, Ak^(A=A) = 0, £k,(A=A) = ^k, ^(A=A) = 0. (21) 

As it turns out, the hybridization Vk is not changed by the renormalization procedure. 



As mentioned before, correlation effects have been neglected in Eq. H16() . First, this means that doubly occupancies 
of /-sites (//^ — //„)(/im — fim) are assumed to be negligibly small though they are not properly excluded by the 
choice of uncorrelated / operators in 7io,A- Note that this assumption is also used for the subsidiary condition within 
the SB approach 7]. Doubly occupied /-sites could in principle be generated by the non-local /-part of Hq.x. As it 
turns out, this is explicitely excluded by keeping the correlations in the interaction part 7ii.\. 

For the following, the commutator of the unperturbed Hamiltonian Ho, a with the operator fl ^Ckm has to be 
evaluated. By introducing the unperturbed Liouville operator Lo,a, which is defined by Lq^a^ = ['Ho,a,^] for any 
operator variable A, one finds 

Lo,A/k„iCkm = [Wo,A,/k„Ckm] = (22) 

= (£/,A - ek,A) fLc^rn + ^ E (^ ~ ^^^)^^'^>^ e'(k-k')R, e*'^'/L^.™Ck™ 

l\f3/2 / , \ "mm';^-»k',A c e J imJ jmJ jm"^\s.m 

where contributions which lead to doubly occupied / sites have been neglected. The second and the third term on 
the r.h.s. of Eq. (|22|) follow from the special form of the anticommutator relations Q. Obviously, only / electron 
operators belonging to different sites i y^ j enter the second term on the r.h.s. of Eq. H22|l . Therefore, as approximation 
one may replace the operator Vjm by its expectation value 

D = (!?,„) = 1-^^I-I(nj) (23) 



"f> - TAfLM (24) 



where 

m 

is the averaged occupation number of / electrons at site j. Note that D is independent of / and m. Furthermore, 
^ fkm i^ replaced by fl.^. Finally, we neglect the third term on the r.h.s of Eq. l(^ . which represents spin- flip 
processes, and all contributions leading to doubly occupied / sites. (Similar approximations will also be used later.) 
Consequently, Eq. (|22|l simplifies to 

I^O,A/k,„Ckm = (e/,A + Ak,A - Aa - ek,A) /k^Ckm (25) 

where 

^^ =- ^E^k,A (26) 

k 

is the averaged / dispersion. Thus, f^j^cum is an approximate eigenvector of the Liouville operator Lq.a- The 
corresponding eigenvalue is the excitation energy Sf^x + Ak.A ~ A^ — £k,A- Furthermore, Eq. H25|l can be used to 
evaluate the action of the projector Pa in (|17|l so that Hi.x can be rewritten as 

Wi,A = PaWi = E ® (^ - 1^/.^ + ^k.A - Aa - Ek.Al) Vk (/Lckm + h.C.) (27) 

k,T?i 

where the ©-function restricts the particle-hole excitations to transition energies smaller than A. 

B. Generator of the unitary transformation 

In the next step let us evaluate a new effective Hamiltonian 7i(A-AA) which is obtained by a further elimination 
of excitations within a small energy shell between (A — AA) and A. According to 114|) 7i(A-AA) is obtained from an 
unitary transformation 

H(A-AA) = e^^-^^^Ae-^^-^^ (28) 



where X\,ax is the generator of the unitary transformation. For the exphcit form of Xx^ax let us make the fohowing 
ansatz 



X 



X.AX 



J2 ^k(A, AA) ek(A, AA) ifl^c^ra - cL/k™) 



(29) 



where 0k(A, AA) is the product of two 8-functions 



ek(A,AA) = e(A-|e/,A + Ak,A-AA-£k,A|) 

X0 [|e/,(A-AA) + Ak,(A-AA) - A(A_AA) 



£k,(A-AA)|-(A-AA)] 



(30) 



The operator form of X\^ax is suggested by its first order expression which is easily obtained by expanding 1)14(1 
in powers of Tii and using Eq. l(TB|) (compare Ref. §). The yet unknown prefactors Ak(A, AA) wiU be specified 
later and depend on A and AA. It wiU turn out that ^k(A, AA) contains contributions to all powers in Vk- Note 
that the ansatz ((28|l also corresponds to the operator structure of the exact generator of the uncorrelated Fano- 
Anderson model (see appendix A). We expect the ansatz (|29|l to be a good approximation also for the correlated 
model in which conduction and localized / electrons strongly couple. Finally, the 8-functions in H29|) result from the 
restriction of TCx to particle-hole excitations with \ef^\ + Ak,A — Aa — Ek.Al < A and from the corresponding restriction 

A(A-AA) ~ £k,(A-AA)l > A — AA for the renormalizcd model Tiix-AX)- The two 0-functions 



^/,(A-AA) 



A- 



k,(A-AA) 



in 6k(A, AA) confine the allowed excitations. 

To determine the unknown parameters Ak(A, AA) of the unitary transformation [compare H29|l ] we will use Eq. H15|l . 
First, we have to carry out the unitary transformation H28|l explicitly 



n 



(A-AA) 



£AAEe''''"'/k™/k -"""^'"^ 



^kme 



k.TTl 



^ £k,A e^^-^^4„Ck™e-^^'^^ 



k,T?i 



Zl^k,Ae 

k,m 

E^k. 
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Xx 



J V.m J k 



'km./km I 6 

' NL 



-Xx 



(31) 



.Xx 



/km'^k 



km^km 



.-Xx,. 



+ Ex. 



The transformations for the various operators are given in appendix El [see Eqs. l|B33|l - l|B36p ]. For instance, the 
transformation of c]^,yiC\im reads, (|B^ 



„Xx,^X„\ „, „— ^A.AA 



'^km'^kJi 



— ek(A,AA)A(A,AA)(/^„Ck„ 



^-D- E UU 



m{^m) 



(32) 



^ ek(A, AA) {cos [2VdAu{X,A\)\ -i}\d (/^„Jk 



NL 



[jmJn 



Dclci,„r - (cl 



'km'^km 



1-^- E (4/..), 



fn{^rn) 



+ -7= Ok (A, A A) sin 2Vi5Ak(A,AA) < /L.Ck™ + h-c 
2\/ D 



^ (/kmCkm + h.' 



l-D - 



fn{^rn 



yJm-l>nJ 



Similar expressions can also be found for the transformation of the remaining operators. In deriving l|B33|l - l|B36|l an 
additional factorization approximation was used in order to keep only operators which are bilinear in the fermionic 
creation and annihilation operators. Spin-flip contributions have been neglected. Moreover, it was assumed that the 
number of k points which are integrated out by the transformation from A to (A — AA) is small compared to the total 
number of k points. 



As already mentioned in the introduction, an expansion with respect to the degeneracy Vf is utiHzed in the slave- 
boson mean-field (SB) theory jaQ- ^^ contrast, here we incorporate all '^/vf corrections which will be reflected by 
the expectation value D as defined in Eq. H23(l . As one can see from the two terms of the anticommutator of Eq. ((SJ, 
new renormalization contributions are included by which a localized electron at an occupied / site is annihilated and 
instead a conduction electron is created. These processes are of order 1/vf smaller than the usual processes included 
in the SB theory by which a conduction electron is annihilated and instead a localized electron is generated at a 
formerly empty / site. 

In the next step let us determine the unknown parameters Ak(A, AA). For that purpose, we insert 7i(\-A\) from 
(|31|l into Eq. 115|) and use the transformed quantities (|B33() - (|B36ll . Thus, one finds 



ek(A,AA)tan[2\/:D^k(A,AA)l = ek(A,AA) ^z! ^\ x • (33) 

L J £f,\ + D (Ak,A - Aa) - ek,A 



The condition 133() for j4k(A, AA) guarantees that 7Y(a-aa) does not contain matrix elements with transition energies 
larger than (A — AA). Obviously, there is a strong similarity between H33I) and the corresponding result (JAIII) of the 
Fano- Anderson model. However, the generator H33|l of the PAM contains some deviations which refiect the influence 
of the strong electronic correlations at / sites. It is important to note that the expression 1^^ for 74k(A, AA) is 
non-perturbatively in Vk and is not restricted to some low order in Vk. Moreover note that the values of Ak{X, AA) 
are determined by Eq. (|32|l only for the case that the excitation energy (e/.a + Ak,A — Aa — ek,A) fits in the energy 
shell restricted by 0k(A, AA). For all other excitations Ak(A, AA) can be set equal to zero. Thus, the parameter 
Ak(A, AA) of the generator ATa^aa is given by 



Ak(A,AA) = <^ 2^ 



arctan 



2VDVu 



Ef^x + D (Ak,A - Aa) - £k,A 



forek(A,AA) = l ^3^^ 



forek(A,AA) ==0 



C. Renormalization equations 

In the following we derive the renormalization equations for the parameters of the Hamiltonian. For that purpose 
we compare two different expressions for 7i{x-A\) ■ The first one is obtained by rewriting the renormalization ansatz 
[Eqs. ^ and ^] at cutoff (A - AA) 



"^(A-AA) - e/,(A-AA)2j'^km/km + 2j^k,(A-AA) (/k,„/krn)^^^ (35) 

k,m k,r?i 
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X! £k,(A-AA) CkmCkm + £^(A-AA) +^VwPx [/k„iCkm + h 
k.,rn k.m 



The second equation for H(x-AX) is found from Eq. (|^ by inserting HB33p - (|B36|l . By comparing in both equations 
the coefficients of the operators cj^^^Ckm, ( /km/km ) , and {fmfni)^ we find the following relations for the parameters 



at cutoff A and (A — AA) 



£k,(A-AA) ~ ek,A 



i [e/,A + D (Ak,A - Aa) - ek,A] {cos [2VdA^{X, AA)] - l} 



- VdVu sin \2VdAi,{X, AA)j , 



(36) 



^k,(A-AA) - ^k,A 



[^k,(A-AA) -Ska] 



1 1 



/,(A-AA) - e/,A = -7^^XlK,(A-AA)-£k,A] 1 + (i^/ - 1) (4 



^/-l 1 



"km'-kn 



(37) 
(38) 



41)3/2 



nY.{ [^/^^ + ^ (^k,A - Aa) - ek,A] sin [2\/:DAk(A, AA) 



- 2^Vk {cos [2V:DAk(A, AA)] - l}} (/Lck™ + h 

^ ^ E [^/-^ - ^ (^k,A - Aa) - £k,A] Ak(A, AA) 



X (/l„ck„ + h.c.), 



l-D vs 



Eix-AX)-Ex = -{l-D)—^N[ef^^x-AX)~£f,x] 



D 



^! 



lYl K,(A-AA) -£K\] 



(39) 



Note tfiat Eq. H39(l for the energy shift Ex foUows from the comparison of the remaining c numbers. The renormahza- 
tion equations J^ - (|^ still depend on the expectation values (cj(.^Ckm), {ftm'^km + h.c. V and D [compare l|^ ] 
which have to be determined first. In the following we will discuss an approximate evaluation for these quantities 
which enables us to solve the renormalization equations (|36ll - (l39ll . The parameter Ak(A, AA) is given by H34|l . 

A factorization approximation was employed above to reduce all renormalization contributions to operator terms 
which are bilinear in the fermionic creation and annihilation operators. In principle, the expectation values are 
defined with the equilibrium distribution of Tix since the renormalization step was done by transforming Tix to 
'^(A-AA)- '^A still contains hybridization terms between / electrons and conduction electrons [compare (| 17(1 ] so that 
there is no straight way to evaluate the expectation values. There are several approximations possible to circumvent 
this difficulty. Firstly, one can evaluate the expectation values by using the unperturbed Hamiltonian Tio.A which is 
diagonal. In this approximation the renormalization equations H36|l - H39|l can easily be evaluated numerically. The 
result for the renormalized / level found in this way is in good agreement with the slave-boson mean field result. 
However, the quasiparticle energies turn out to be discontinuous as function of k. This behavior has to be interpreted 

as an artifact of this simple approximation and is caused by the vanishing of expectation values (/kmCkm + h.c.\. 
The second possible approximation is more difficult and consists in calculating the expectation values with respect 
to the full Hamiltonian H instead of Hx- In this case the renormalization equations (|36|l - (|39|l can not be evaluated 
directly because the expectation values are not known. The starting point to find these quantities is the relation 



{A) = 



Tr{Ae-'^'^) _ Tr{Axe- 
Tr (e-/3«) 



-/3Ha 



Tr (e" 



-pn^ 



(40) 



which follows from unitarity (for any operator A). By setting up additional renormalization equations for the trans- 
formed operators Ax one can determine the expectation values (AY Note that in the equations for ^a the unknown 
expectation values enter again so that they have to be solved self-consistently. This approach is rather involved but 
has the advantage that expectation values in ((36(1 - ((39|l no longer depend on the cutoff energy A. 



Renormalization equations for transformed operators also have to be used if dynamical correlation functions are 



evaluated. For example, to find the densities of states of the / electrons 

-"/M = J^J2([fL^^i^ + ^)hrn\) (41) 

km \ ' 

and of the c electrons 

= (^) = ^E(^L>'5(L + ^)ckm]^) (42) 



P< 

iV 

kr?i 

one has to apply the renornialization transformation on f^^ and cj^.^. This will be done in appendix O Note that in 
Eqs. H41I) and (|42|l the Liouville operator L of the full Hamiltonian was introduced which is defined by hA ~ [HjA] 
for any operator variable A. 

IV. ANALYTICAL SOLUTION 

Alternatively, one can also find approximate analytical solutions for the renornialization equations H36(l - (|39|l . For 
this purpose three approximations have to be used: 

(i) All excitation values are calculated using the full Hamiltonian H (see the discussion at the end of Sec. llllll . i.e. 

(■••>«.«(•■■)« = (••■>, (43) 

so that they are independent from the renornialization parameter A. 
(ii) The A dependence of the renormalized / level will be neglected, 

^/,A « £/■ (44) 

The spirit of this approximation is similar to that used in the slave-boson theory. There a renormalized / 
level is used from the very beginning. Within the present treatment one might expect that Sf,\ increases with 
decreasing A from its initial value Sf and reaches its final value if already at finite A. 

(iii) The averaged dispersion of / electrons will be neglected, 

k 

These approximation enable us to map the renornialization equations of the PAM to those of the exactly solvable 
Fano- Anderson model (see appendix fX|l . 

A. Quasi-particle energies 

Eq. (|37|l depends on differences of the parameters of the transformed Hamiltonians at A and (A — AA). Therefore, 
this equation can easily be integrated between a lower cutoff A ^ and the cutoff A of the original model. One finds 

Ak = - [£k - £k] . (46) 

where the initial parameter values (|21|l were used. Furthermore, we have defined A^ = A^ rA^o) ^-i^d et — £k,(A^o)- 
Eq. H39|l can also be integrated in the same way so that we find 

E ^ - {n{)N {if ^ef)-^Y. ^^^ ' ^k) • (47) 



k 



Here, again the initial parameter values (|21l) and Eq. (|23|l have been used. Furthermore, we have defined E = E/^x^o) 
and if — e^^ (a_+o)- The second term on the r.h.s. of 147(1 vanishes if we use Eq. (|46|l and approximation (|45() so that 
we obtain 

E = -{n{)N{if-ef). (48) 
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It is more difficult to solve the remaining renormalization equations (|36|l and H38|) . First, the approximations H43|) - 
(|45|l lead to a decoupling of Eq. H36|l for different k values. Thus, by eliminating excitations from large to small 
cutoff values A each k state is renormalized only once. Such a step-like renormalization behavior is also obtained 
in the Fano- Anderson model (see appendix fX|) . A further similarity to this simpler model is found by inserting the 
approximations gSJ-gSj into Eqs. ^^ and ^^ 



i.n(2^AA = I^^ 



£/ 



1 



(49) 



e^-e^ = 2 



cos h^DA^ - 1 [£k - £/] - VDVk sin hVOAy,) . (50) 

Here, the step-like renormalization behavior and the initial parameter values H21|l have been used. Note that in the 
case of a step-like renormalization the parameters ^k of the generator of the unitary transformation do not depend 
on A and AA. The above equations are very similar to those obtained for the Fano-Anderson model [compare with 
Eqs. (jAlip and ljA10|) ]. In particular, the equivalence of the one-particle energies can be seen by replacing \fDV\^ 
by 14c- Moreover, one immediately finds from l|36(l and (|5()|l the following result for the renormalized c electron 
one-particle energy 

~e. ^ ^ ~ '^h_lAw, (51) 

where 



Wk = ^(£k-£7)'+4i?|ykP. (52) 

Obviously, Eqs. (|46|l and H51() also determine the / type quasi-particle excitation energy which is given by 

r, • ? ^A £/ + £k , sgn(£/-£k) ,„ 

^k .= £/ + Ak = 1 Wk (53) 

where approximation (|45|) has been used. Thus, we have obtained two quasi-particle excitations. According to Hl()|) 
and H53|l the renormalized Hamiltonian reads 

n := 7i(A^0) = E '^k /il„jkm + Y. ^~k cLckm + E. (54) 

k,rn k,m 

Note that in 7i the hybridization was completely used up for the renormalization of the parameters of Tip. Also the 
eigenmodes /km and Ckm of 7i do not change their character as function of the wave vector due to the presence of the 
sgn- functions in H51(l and (|53|l . Instead they remain /-like or c-like for all values of k. Furthermore, the one-particle 
energies ()51|l and (|53|l still depend on two unknown quantities, namely, the renormalized / level £/ and the expectation 
value D [see Eq. (|23|l ] which have to be determined in the following. 

B. Free energy and equations of self-consistency 

First, let us calculate the averaged / electron occupation number (nf) from the free energy F. Note that 'Hx is 
connected with the original Hamiltonian 7i by an unitary transformation. Therefore, the free energy can also be 
evaluated from l-i\. In particular, the relation 

F = --InTre"'^^ = -- InTr e"'^^ =: F (55) 

P P 

holds. Because of TL describes an non-interacting Fermi system the free energy F can be easily calculated 

F = -^ ^ In [1 + e-^'-] - ^ 5] In [1 + er^'^-] + E. (56) 

^ k ^ k 

The / electron occupation number is found from F by functional derivative 
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which can be easily performed. We finally obtain a relation of the following structure 





{...}(§)+{...} 



'd{h{) 



de 



f 



(58) 



We are interested in solutions of the rcnormalization equations which describe mixed valence and heavy Fermion 
behavior. For these cases the derivatives in Eq. (|58|) are non-zero so that solutions can be found by setting both brace 
expressions equal to zero. In this way the following self-consistent equations for the renormalized / level if and the 
averaged/ electron occupation number {n{) are found, 



iE/(^~0{^sgn(.,-e.)£|^} 



'iJlfi^^^U+'Snie^-^ff-^ 



(59) 



^/-£/ = 



N 



J2 ^g"( 



£/ 



c)/(ek) 



l^kl^ 



N 



Esgn(ek-£/)/(wk) 



Wk 



(60) 



Note that these equations are quite similar to those which are found in the slave-boson (SB) formalism Q- In 
particular, the limit i// — > oo of Eqs. (|59|l and l|6()|l perfectly agrees with the SB equations. 



C. Expectation values 

The remaining expectation values (cl^„^Ckm) and (/j[„jCkm + h-C-) can also be evaluated from the free energy (|55|) 



^km'-^kn 



1 dF _ 1 &H 
vf (9ek vf \ dew 



H 



(61) 



/Lckm + h.C. 



1 dF 



"f 



dVi. 



1 dn 



H 



Both expressions can be evaluated similarly to H57|) . By using H59|l and (|60|l we find 



'^km'-kr] 



1 - Sgn(£/ - £k! 



Ek-e/ 



1 - sgn(ek - if) 



2VFk 

£k — £/ 
2Wk 



/(£k) 
/('^k), 



(62) 



(63) 



/kmCk™ + h.c. ) = - 2 sgn(£/ - Ek) -^ fiik) - 2 sgn(£k - £/) -^ fi^^k)- 



Note that also Eqs. (|63|l and H64() are very similar to the corresponding SB results. 



(64) 



D. One-particle operators and density of states 

Next we calculate the densities of states of the / and c electrons [compare Eqs. (|41|) and (|42|l ]. For that purpose, we 
have to integrate the rcnormalization equations ljC9p and ICIOII to determine the transformed one-particle operators. 
In the case of the analytical solution, ljC9|) and (|C10|) can be exactly solved if the basic approximations H44|) and 145|) 
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are used. As already discussed above, in this case the different k values are not coupled with each other and we obtain 
a step-like renormalization behavior. Thus, we find 

Uk — cos(vI)Ak) and Wk = ^^sin(vZ?^k) (65) 



D 

where the initial parameter values (|C2(I were used. Furthermore, we have defined Wk = UkAx^o) ^-nd Uk — I'kXA^o)- 
Combining the generator of the unitary transformation (|49|l , the normalization condition (|C3|) , and Eq. H()5|) we finally 
obtain 

\u^.f = ^{l-^l^«gn(^~/-^k)} (66) 

l^.|^ = ^{l + ^sgn(s,-Sk)}. (67) 

Thus, the coefficients -Uk and Wk can be directly calculated from the results of the self-consistent equations H59|l and 

To calculate the densities of states (|41f) and (I42|l we use the relation 14U|) which follows from the unitary of all 
operators. Thus, Eqs. (|41|) and 142|) can be rewritten as 

Pf 



(^) = ^ E ( [fLM ^ 0),5 (l + cj) /k„.(A ^ 0)] ^\ , (68) 

km ^ ' ^ 

(^) = ^ E ( [irn ( A ^ 0) , <5 (l + c.) ck™ ( A ^ 0)] \ . (69) 



Pc 

where L is the Liouville operator of the final Hamiltonian H which is defined by hA = [HjA] for all operator variables 
A. Due to the structure (|54|l of the final Hamiltonian Ti, Eqs. Ht)8|) and (|69|l can be easily evaluated. Using ljCl|l . 
(EH), lEH), and (EH) we obtain 



(io) = I?^^{|uk|'(5(c^-c^k) + C|«k|'<^(c^-ek)}, (70) 

H = ^E{l^k|''5(^-ek) + I?|i)k|'<5(c^-c:;k)}. (71) 



E. Results and comparison with slave-boson mean-field theory 

In this subsection we shall compare the results of our analytical solution discussed above with those of the slave- 
boson mean-field (SB) treatment. As already mentioned, the limit i/^ — > oo of the derived self-consistent equations 
(|59|l and HtJUfl is completely equivalent to the SB equations. Furthermore, in this limit the expectation values of 
our analytical solution [see Eqs. Ht)3|) and Ht)4|) ] and the SB treatment also perfectly agree. Therefore, we want to 
concentrate on the case of small degeneracy Vf . At this point it is important to notify that we have never exploited 
an 1/i^f expansion in the derivation of the analytical solution of the PAM so that it is valid for large as well as small 
degeneracy Vf. 

For simplicity, let us consider an one-dimensional PAM with 50000 lattice sites, a linear dispersion relation for the 
conduction electrons, and a k independent hybridization V^ — V and compare our results with those of the slave-boson 
mean- field (SB) theory. In particular, we are interested in the dependence of the results on the degeneracy Vf. 

At first, let us fix the degeneracy of the angular momentum to Vf = 4. The other parameters are chosen as follows 
i^fV^ = 0.36, Ef = —0.3, chemical potential /x = 0, and T = 0.00001 where all energies are given in units of the half 
bandwidth. As can be seen from Fig.^ the renormalized quasi-particle energies, i.e. Sk and ujk, obtained from l|51|l . 
(|53|1 . H59|l . and H60|) (full and dashed thick lines), and the quasi-particle bands of the SB theory (dotted lines) seem 
to be quite similar. However, the averaged / occupation (nf) = 0.855 and the renormalized / level if — 0.071 differ 
significantly from the SB results ((n^) = 0.586 and if = 0.115). These differences are mainly caused by the fact that 
we have taken into account all l/i^/ corrections which are absent in the SB treatment. Note that l/vf corrections 
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FIG. 1: Dispersion relations of an one-dimensional PAM (iV = 50000, Vf = 4, VjV'^ = 0.36, e/ = -0.3, ^ = 0, T=0.00001). 
Here, the unrenormalized one-particle energies Sk and £/ are plotted with full and dashed thin lines. The renormalized 
quasiparticle energies ek and iDk are shown with full and dashed thick lines. Furthermore, the quasiparticle energies of the SB 
approach are drawn by use of dotted lines. 
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FIG. 2: Density of states of the / electrons (upper panel) and of the c electrons (lower panel) where all parameters are chosen 
as in Fig. Q A broadening of the (5-functions of 0.0001 is used. The results of the analytical PRM solution (SB theory) are 
drawn as solid (dashed) lines. 
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FIG. 3: Dependence of the renormalized / level if (left panel) and of the averaged / occupation {n{) (right panel) on the 
degeneracy Vf of the angular momentum where all other parameters are chosen as in Fig. Q The results of the analytical 
solution (SB theory) are drawn using solid (dashed) lines. 
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FIG. 4: Dependence of the renormalized / level if (left panel) and of the averaged / occupation (n{) (right panel) on the 
original / energy Sf where all other parameters are chosen as in Fig. The results of the analytical solution (SB theory) are 
drawn using solid (dashed) lines. 
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FIG. 5: Averaged / occupation (n[) as a function of temperature T where all other parameters are chosen as in Fig. The 
result of the analytical solution (SB theory) is drawn using solid (dashed) line. 



allow for additional renormalization processes which lead to a lowering of the free energy for the whole parameter 
space. 

It is well known that the quasi-particles of the SB theory change their character as function of k between a 
more /-hke and more c-like behavior. As was mentioned before, in the present treatment excitations do not change 
their character as function of k. However, the quasi-particle energies show jumps in their k dependence where the 
renormalization contributions change their sign from positive to negative values or vice versa. Note, that the various 
parts of the quasiparticle bands fit perfectly together (see Fig. P). 

As compared to the dispersion relations plotted in Fig. 2] the densities of states of the /- and c-electrons in Fig. |21 
show much better the differences between the results of the present analytical solution and of the SB treatment. In 
particular, the smaller value for the renormalized / level if obtained from our PRM treatment leads to much higher 
density of states at the Fermi surface than the SB treatment. Note that such an enhanced density of states at the 
Fermi energy is a clear signature of heavy fermion behavior. 

Next, we discuss the dependence on the degeneracy parameter Vf. For that purpose we vary Vf by keeping 

VfV'^ = 0.36 fixed. In contrast to the SB results for (nf) and if, which are almost unchanged (see Fig. O, the 
analytical solutions show a remarkable dependence on the degeneracy Vf. In particular, for small Vf, the l/vf 
corrections included in the PRM approach lead to serious deviations from the SB results. From these additional 1/vf 
corrections follows a more pronounced heavy Fermion behavior. As already mentioned above, the limit z/j — > cxd of our 
analytical solution perfectly agrees with the SB theory. To perform this limit one has to replace the expectation value 



i^by(l 



^/ 



)) so that all processes are neglected by which a localized electron at an occupied / site is annihilated 



and instead a conduction electron is created [compare discussion below Eq. (|32|l ] 

In Fig. 21 the renormalized / level if and the averaged / occupation (n| ) are plotted as functions of the original / 
energy £f. The momentum degeneracy has been fixed at Vf = 4. As is seen, the l/vf corrections do not only cause a 
dependence of the results on the degeneracy Vf (as shown in Fig. (JJ but also lead to a reduction of the stability range 
of heavy Fermion type solutions. 

The \/vf corrections also affect the thermodynamic properties of the system. In Fig.|Slthe temperature dependence 



of the averaged / occupation (n^) is shown where Vf has been fixed to Vf = 4. We observe that 



goes with 



increasing temperature much faster to 1 than the SB results. Thus, the ^/vf corrections lead to a lowering of the 
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Kondo temperature Tk which may be defined as that temperature at which (n^) becomes 1. 

As can be seen from Fig. 0] the analytical PRM solution breaks down when the unrenormalized / level becomes 
smaller than some critical values. A similar behavior is also known from the SB solution. For instance, the solution 
for {ni) breaks down when for fixed £f the chemical potential fj, is increased beyond some critical value |l2l |. The 
reason for this breakdown is not completely clear. May be, it is due to some rough approximation used both in 
the PRM treatment and the SB theory, for instance those from Sect. III. A which have their counterparts in the SB 
treatment. Alternative, the breakdown of the PRM and the SB solutions might be a signature of a genuine phase 
transition. Recently it was suggested |l3ij that in certain systems like CeCu2Si2 there might be a transition between 
an intermediate valence regime with fluctuating / charges and a regime with integral / charge when the pressure is 
decreased. In the integral regime which is described by the Kondo Haniiltonian there is no longer a renormalized 
/ level at the Fermi level. This might be the reason that the self-consistent solution for if no longer exists. Note 
however that such a phase transition does not appear in a recent alternative discussion of the PAM on the basis 
of Hubbard operators in Ref. 12?. This approach is based on an extended chain approximation and gives the same 
quasi-particle energies (|51|l and H53|) . However, it leads to completely different equations for the renormalized if level 
and for the averaged / occupation {n{ ) . Results have been found which are very similar to the SB solutions for those 
parameter regimes where the SB solution exists. In contrast, the PRM solution leads to substantial deviations from 
the SB results in particular for small values of i//. Note however that apart from the l/i^f corrections we have used 
similar approximations as in the SB theory to derive our analytical solution. 

V. CONCLUSIONS 

In summary, in this paper we have applied a recently developed projector-based renormalization method (PRM) to 
the periodic Anderson model (PAM) in the limit of infinitely large Coulomb repulsion at / sites. By using an additional 
factorization approximation we have derived renormalization equations for the parameters of the Haniiltonian. In this 
way, the PAM is mapped to a free system consisting of two uncorrclated quasiparticle bands. Similar uncorrelated 
Hamiltonians have been also derived before by different theoretical approaches, such as the Gutzwiller projection _a| 
and the SB theory [y, |3 where l/i^f expansions have been exploited. In contrast, the present approach is valid for 
any Vf. Due to the factorization approximation certain expectation values enter which prevent a direct numerical 
evaluation of the renormalization equations. In principle, the expectation values could be determined self-consistently 
by deriving additional renormalization equations also for the operators which enter the expectation values. This has 
not been done in this paper. 

To obtain instead an analytical solution we have used a renormalized / level which was assumed to be constant 
during the renormalization process. The spirit of this approximation is similar to that used in the SB theory [y, 0] ■ 
We obtain self-consistent equations for the renormalized / level if and the averaged / occupation (n{) which are 
quite similar to those of the SB theory [3. In particular, in the limit ff —^ 00 our solution perfectly agrees with the 
SB result but strongly differs from it for smaller values of Vf. To compare our results in more detail with those of 
the SB approach we have also considered an one-dimensional PAM with a linear dispersion relation of the conduction 
electrons and a k independent hybridization. Note that the character of the obtained two quasi-particle bands of 
the two treatments differ. Whereas the quasi-particles of the SB theory change their character as function of k from 
a more /-like to a more c-like behavior and vice versa the excitations of the PRM treatment do not change their 
character. Instead, the quasi-particle energies show jumps as function of k where, however, the various parts of the 
quasi-particle bands fit perfectly together. The influence of the degeneracy i/f has been studied by varying i/f with 
fixed VfV^. Whereas the SB results are almost unchanged, our analytical results show a remarkable dependence on 
the degeneracy Vf. Especially, serious deviations are found for small values of Vf. 

Finally, from a more technical point of view, note that in in this paper the PRM method was applied to a physical 
system for the first time without using any perturbation theory. 
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APPENDIX A: RENORMALIZATION OF THE EXACT SOLVABLE FANO-ANDERSON MODEL 

In this appendix we illustrate the usefulness of the projector-based renormalization method (PRM) for the case of 
a simple model. We apply the approach of Sec.mto the exactly solvable Fano- Anderson model [3,I3- This model was 
already discussed in the framework of the present approach in Ref. |3- However, now this will be done in a consequent 
non-perturbative manner. 

1. Model 

The Fano- Anderson model consists of dispersionless / electrons which interact with conducting electrons. Thereby 
all correlation effects are neglected. The Hamiltonian reads 

n = Ho+Hi, (Al) 

Wo = ^ (e/ /km/km + Ek C^^ra^km) , 
k,m 



Hi = V Mc (/Lckm + 4„/km 



As in Eq. (^ the index i denotes the / sites, k is the wave vector, and Vk describes the hybridization between 
conduction and localized electrons. The excitation energies £k and e/ for conduction and localized electrons are 
measured from the chemical potential fi. Both types of electrons are assumed to have the same angular momentum 
index m with values to = 1 . . . ^'/. Of course, the model is easily solved and leads to two hybridized bands 

^ = E 4"^«L„«k™ + E ^i^^f^LP^rn (A2) 

k,T7i k,m 

where 

(a,/3) _ £k + £/ 1 

^k - -^±^Wu, 



\2 , ^ITA |2 



The eigenmodes a]^^ and /3^„j are given by linear combinations of cl^^^^ and /i[„j, 

^km — ""k/km + ^kCjj.„j, p,j.„j = — ^k ./km + "k Cj^.,. 



ai 



12 _ Ml £k-£/\ , , ,2 _ 1 /, £k-£/ 



"kl^ - 7:^--^TTr^ U M' = 7T 1 



(A3) 



2 V Wk / 2 V VFk 

2. Renormalization ansatz 

In the renormalization approach we integrate out particle-hole excitations of conduction and / electrons which enter 
due to the hybridization term Tii. We expect to finally obtain from the PRM an effectively free model. The starting 
point of the method is a renormalized Hamiltonian Tix which is obtained after all excitations with energies larger 
than a given cutoff A have been eliminated. Due to the result of the preceeding section it should have the following 
form 

-Hx = Ho^x+Hi^x, (A4) 

Ho, A = 2^ (^£k,A fkmf^m + £k,A ^kmCkm j , 
k,?n 

Hi,A = E ^k (/Lck™ + cL/k™) = Pa Hi . 



k,771 
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As it turns out, no renormalization of the hybridization Vk occurs so that Vk is assumed to be A independent from 
the beginning. Like in the exact diagonahzation different k states do not coupled during the renormahzation process. 
Thus, by ehminating excitations from large to low A values, each k state is renormalized only once leading to a step 
like renormalization behavior. This means, for a given cutoff A all k states with excitations |ek — £/| > A have already 
been renormalized whereas those with jet — £/| < A have not. Thus, TYa can be written as a sum of two parts 
nx=n^ + Til where 



n = 


k,m 
|^k-^/|<A 


-nl = 


k.m 
|£k-e/|>A 



^/ JkmJ^"i + ^k C^^rn'^km + 14; ( /km'^km + Ckm/kmJ 



[fk JkrnJi^ 



k '^km'^kr) 



(A5) 



(A6) 



7i^ is the unchanged part of Tix whereas Ti^ is renormalized due to the elimination of excitations |ek — e/| larger 



than A. ej^ and e^ denote the renormalized energies. 



Transformation of the Hamiltonian 



For the explicit evaluation of 7i^ let us apply the unitary transformation 1)11(1 on the original Hamiltonian Ti. 

Ux = e^^Ue-^K (A7) 

For the generator Xx of the unitary transformation an exact expression can be given. By inspection of the perturbation 
expansion in terms of 14; [8| one finds that Xx must have the following operator structure 



Xx = 



E ^k {fL 



.^kn 



^kmJ^ 



lz.ni 



(A8) 



k,rM 

|ek-£/|>A 



with yet unknown prefactors Ai^. Eq. (|A8|) will be taken as ansatz. Then Ti^ can be easily evaluated since only k 
values with jet — £/| > A renormalize the Hamiltonian. Due to the fermionic anticommutator relations different k 
values are not coupled. To find 7i^, we consider the transformation of the various operators. For instance, we obtain 



^ '^km'-km 6 Cj.„jCkm 



(A9) 



= ^0 (kk - e/l - A) {[cos (2Ak) - 1] (c^Ck™ - fLh 

+ sin (2Ak) (/k^Ckm + 4„/kmj I . 
Similar relations can also be found for the transformations of the operators /k,„/km and (/kmCkm + cj(.„/km). Thus, 



H^ reads 



n>x 



E 



k,m 

ek-e/l>A 



£/ - 2 [cos(2v4k) - 1] (£k - £/) + ^k sin(2Ak) 



£k + 2 [cos(2Ak) - 1] (£k - £/) - Vk sin(2Ak) 



/km/k 



Ckm^kn 



(AlO) 



Vk + 2 sin(2Ak) (£k - £/) + ■^i [cos(2A) - 1] 






^U.n 



1 ,7krn ) 



In contrast to the expected form (|A6() for 7i^ the expression HA10() still contains a hybridization part proportional 

to (/i[„jCkm + Ck„j/km) with excitatiou energies larger than A. The requirement ({T^. QaHa — 0, leads to the following 
condition for A^ 

2Vk 



tan(2^k) 



Ef -£k 



(All) 
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Eq. (|A11|) guarantees that the hybridization vanishes in ljA10|l and TC^ becomes diagonal. Note that according to 
(|A11() , the quantity A]^ changes its sign when the energy difference Sf — Ek changes its sign. By inserting (|A11|I into 
l|A10|l one finds 

^A = E (4/L/k™ + £lcLck™) (A12) 

k.m 

£k-£/|>A 

where the renormalized energies are given by 

-/ _ £/ + £k , sgn(£/-£k) 



2 2 



Wk, (A13) 



.c _ £/+£k sgn(£/-£k) „^ 
£k - -^ ^ Wk. 

For A — > the Hamiltonian is completely renormalized. The final Hamiltonian Ti. :— 7i(^^o) reads 

^ = E (^t /La™ + Ek cLck™) . (A14) 

k,?n 

Note that the final result (|A14|I corresponds to the diagonal Hamiltonian of ea. (|A2|l . In particular, all expectation 
values completely agree between the two approaches. (To calculate expectation values by using Ea. ljA14|) one also has 
to transform the operators which enter the expectation values. For more details see Ref. |^.) However, in contrast to 
the eigenmodes a]^^ and /3^^ of l|A2|l the present eigenmodes fl^ and c].^ do not change their character as function 
of the wave vector. Remember, ak„ was a more /-like excitation for £k < £/ and a more c-like excitation for £k > £/, 
and vice versa for /3^,„. In contrast, the operators /j^^j^ and Cy.^^ of (jA14p remain /-like and c-like for all values of k. In 
return, the A dependent excitation energies £j^ ^ and £k > show as function of A a step-like behavior at A = |£k — £/|. 
This step-like change guarantees that deviations from the unrenormalized energies £/ and £k stay relatively small for 
all k values. 



APPENDIX B: TRANSFORMATION OF THE OPERATORS 

In this appendix we evaluate the transformation from A to (A — AA) for the various operator quantities of Eq. I|31|l . 
For example 

OO -J 

e^'^'^cLck^e-^^.- = e^^- (cL^k™) = E ^^ ^a.aa (cL^k™) . (Bl) 

n=0 

Here, a new super-operator Xa,aa was introduced which is defined by the commutator of the generator AT^aa with 
operators A on which X^.aa is applied, Xa,aa-4 ~ [X\^^\,A]. Furthermore, let us define a new operator Akm by 

^km = /km/'km ~ Ckm/km 

which is an ingredient of the generator of the unitary transformation 

^A,AA = E^k(A,AA)ek(A,AA)Ak„,. (B2) 

\<.m 
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We have to evaluate various commutators 



"^k'm'j Cj^m'^k 



X, 



k'm' 



\jmJrnj 



= Sk'M Sm'.rn f/kmCkm + h.C 
S,n' 



N 



r" (/k'mCk'm + h.C. ) , 
[Jmlm, J 



X, 



X, 



k'm' 



\^./km/kmJ 
\^/kmJk»nJ 



X 



k'm 



' l^/km-'km 



+ h.c. 



NL 



NL 



NL 



7V3/2 



E {e'''^'e'('''-'^)^^7t„P,„.Ck'™ + h.c.}, 



ij(#i) 



^"'''"' V { J^^^ A^' -'^)^o f^ T) c^' +hcl 

/Y3/2 A^ \^ ^ yim-^jmf-k'm ^^ U.C.J , 






'km-Zkire 



NL 



X 



k'm/ 



'NL 



(B3) 
(B4) 
(B5) 
(B6) 

(B7) 
(B8) 



X, 



, (/k„>Ckm + h.C. 
Xk'm', f/kmCkm + h.C 



^tn' .■ 



— Orri.' ' 



|2^k',k/k„/km - cLnCk'm2?m(k' - k) + h.C. | , 

|25k',k (/km/km + h.C.j - 4,„Ck'mX',„(k' - k) + h.C. | 



(B9) 



(BIO) 



where we neglect all spin-flip contributions. In Eqs. IJB9|I and ljB10|l Fourier transformed quantities are introduced 



A„(k) = l^e^'-^'I?,, 



(Bll) 



Furthermore, we have defined 






L ' X 



NL X 



(^/km/km 



h.C 



NL 



AT Z^'' km /km — j^r /_^ J imJim: 
k i 

^ Z^ JirnJjmf^ — JkmJ^ 

^ E [/L/.™e^'^(^-^^)+h.< 



JmJm I , 



(B12) 
(B13) 
(B14) 



iJii^i) 



We are interested in contributions which renormalize the parameters of the Hamiltonian Tix according to Eqs. pt)|) 
and p7|l . Therefore, an additional factorization has to be carried out in Eqs. (jB6p . IJB7|) . ljB9p . and ljB10|l and to 
keep only operator terms which appear also in Tix. By neglecting more complex operators, namely spin- flip terms. 



21 



Eqs. (|B6l), ||B3, 1191), and (|BTO|) can be replaced by 



-^k'm', (^/km/kmj 



NL. 



= —Sm.m' S] 



'k,k' 



^){-(/i. 



Ck' 



+ (/k'mCk',„+h.C. 



1- D- 



rh{^m) 



ir?ihi 



X, 



k'm 



' \^^km-/kmj 



NL. 



'-^rn.m' I ^k.k' 



^)|^(/k'™Ck'™+h.C.) 



+ (/k'mCk'm+h.C. 



l-Z? 



E 



ifniri 



(B15) 



(B16) 



^k'm', (^/k„Ckm+h.C.j = 2(5„j,,„/(5k,k' < (^/k„j/kmj + {^fLfmj 



Dcl^Ckr, 



'^km'^kri 



l-Z? 



E 



[jmJ'nj 



(B17) 



-'^k'm', ( /kmCkm + hj 



= 2(5m,m'5k,k' S ^ ( /m/m + h.C. 



[jmJmj 



^^km'^km ( ^kmCkm 



l-D- 



7n(7^Tn) 



I 7rh/m I 



(B18) 



where ( /k^'^km + h.c. ) ~ ( /km'^km + h-C. ) has been used. Due to this factorization, certain expectation values 

enter Eqs. (|B15|) - ljB10|) which have to be evaluated separately (compare the discussion in subsection llllCp . By using 
Eqs. (|B3p - IIB5ll and (|B15|) - (|B18|) one finds from ljB2|) for the corresponding commutators formed with the generator 



-'^A,AA,CkmCk 
-'i^A.AA, f/m/r 



ek(A,AA)Ak(A,AA) (/Lckm + h.i 

- ^ E ©k' ( A, A A) Ak ( A, AA) [fl^c^'m + h.c 



[-'^A.AA, {flifm)j] - 



-'^A.AA, ( fmfm 



(B19) 
(B20) 

(B21) 
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X\^A\, (^ /km/km j 



NL. 



ek(A, AA) Ak(A, AA) J D (/^ck™ + h.c. 



+ (/Lckm+h.C. 



1-^- E {flf^), 



m{^m) 



lX]ek'(A,AA)Ak-(A,AA)Ji5(/^,„ 



Ck'm + h.c. 



/k'„Ck'™ + h.c. 



l-D 



rn{^rn) 



Jrn Jrh 



(B22) 



^A,AA, (/L/krn) ^^] = ^ 0k(A, AA) Ak(A, AA) J D [fl^Ckrn + h. 



(B23) 



+ ( /kmCkm + h.c. 



l-D- J2 (/rU 



+ ^ E Qk'(A, AA) Ak' (A, AA) J i? [.fl,„,cu'nr + h.c. 



/k'™Ck',„i + h.c. 



1-^- E fe/- 



7n(^m) 



X 



A,AA, 



(^Jkm/k 



+ h.c. 



NL. 



X 



A,AA, 



(/km/km jj^^ + -'^A,AA, (^ /km/km j 



NL. 



(B24) 



X^,A^JLc^.m + h.c] = 2 ek(A, AA) Ak(A, AA) J {fljlm)^^ + {fUm) 



(B25) 



£'ci„,Ckm- (4 



'km 



Ckrj 



1-^- E (/U-"), 



i(5^m) 



-'^A.AAj/k^Ckm + h.c. 



2 ek(A, AA) Ak(A, AA) { | (/L/k™ + h.c.)^^ + ( /^/„ 



(B26) 



^'^km'^km ( Ci^ Ckm ) 



1-^- E (4/ 



m{^m) 



Therefore, all operators terms appearing on the r.h. sides of Eqs. HB19|) - (|B26|) are traced back to a bilinear form. 
This property will enable us to evaluate higher order commutators with X\^a\ and also transformations like (|B1|I . 
Moreover, we assume that the number of k points which are integrated out by use of the unitary transformation H28() 
is small compared to the total number of k points. This assumption is needed for the evaluation of higher order 
commutators. For instance, the commutators which arise from repeated application of X\^a\ to (/j^^Ckm + h.c.) and 
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(/kmCkm + h-c.) read (n = 1, 2,3, • • • ) 



X^"AA(/Lck™ + h.c, 



= (-l)"ek(A,AA)f2V:DAk(A,AA)l "\(fLc^m + h.c 



1 
2b 



7^ ( /Lckm + h.C. 






^+l /'ft 
,AA l^^k 



A,AA I /km'^km 



+ h.C. 



(-1)" 



ek(A,AA) 2Vi:'Ak(A,AA) 



^Ci^-m.Ckm \Ci,™Ck 



2n+l 



' krr?, -^ ^ 



NL 



JmJn 



-'km'- 



l-D 



E 



JniJrn 



(B27) 



(B28) 



X^"aa ( /Lck™ + h.c. 



7-2n 
''A.AA 



AA \^.'k?i 



,Ckm + h-C. 



Xa"aa (/Lckm + h.c. j + (-1)" ek(A, AA) [^Ak{X, AX) 

X { (/Lckm + h.C.) - (/^„Ck™ + h.c.) } 

(-1) 



2n 



(B29) 



■■A.AA Wk: 



X'"A+M/Lckm+ll.C 



£1 



ek(A,AA) Vi:'Ak(A,AA) 



2n+l 



{(/L/k™ + h.c.)^^-2 



/L/km)^^} 



(B30) 



To trace back all contributions to terms appearing in the unperturbed Hamiltonian 7io,Aj one has to replace all 
Hubbard operators by appropriate expressions in terms of usual Fermi operators. Thus, further approximations 

are needed for the local ( /i[„j/km ) and for the non-local / electron particle-hole excitations I /j[,„/km ) • As 

\ / Li V / JN ij 

was discussed before, due to the strong local Coulomb interaction, only empty and singly occupied / sites are of 
physical relevance. Thus the operator I /j[^/km ) applied on physical states can not generate doubly occupied / 
sites. Therefore, the operator can be replaced by 



y J km J km j ~ \^/km /km j 



(B31) 



The second operator, (./km /km) represents an / electron hopping between different sites without creating doubly 



NL 



occupied / sites. Thus, we may approximate 

(/L/km)^^ - ^ E /i/.™^^''^"--"^' = ^ E ^""/L/.^e^"'"' 

« D (/L/km) 



-R-.) 



(B32) 



i,Ji¥=i) 



'NL 



where in the last equation the creation of doubly occupied sites is only fulfilled within a factorization approximation. 



Finally, by inserting (|B27|) - ljB28|) into IjBip . and by using the approximations (|B31|) and l|B32p one finds 
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^ ''km 'knit; ''Icm'-kni — 



(B33) 



— ek(A,AA)v4k(A,AA) (/^,„ckm + h.c. 






^ ek(A, AA) {cos [2V:DA(A, AA)] -i]\d [flj^. 



NL 



+ \JmJn 



-'^ "-km "^km \ '^km "^km 






1 



2V£' 



ek(A,AA) sin[2V:DAk(A,AA)l <^ f/Lck™ + h-c 



^(/L'^km + li.c. 



1-^- E (4j^ 



Similar equations can be derived for the transformations of the remaining operators 



^Xx 



(^/km/kmj^^ e (^/km/kmj 



NL 



(B34) 



D { e^^'^^c^ckme-^^-^^ - 4^Ckm + ek(A, AA) Ak(A, AA) Ul^c^,^ + h.c. 



1-^- E (/i/'^ 



D 

N 



?|i 



'^ i-k'm^-k'mc ^k'r7i'-k': 



- Gk' (A, AA) Ak' (A, AA) ( /i[,„Ck'™ + h.c. 



1 -D 






yJrhJrfij 



-,X. 



[J ml" 



p-Xx,AX ^ ( f^ f 

, t- \JmJrr 



^n 



^Xx.AX ^'\ p, , ^ — Xx^AX 

^ '-k'm'-k'm'^ 



''k'm'^k'; 



(B35) 



„x. 






Ck'm +h.c. e 



-X. 






Ck'm +I1.C, 



ek(A, AA) {cos [2VdAi,{X, AA)] - l} J (/Lck™ + h.c 



-R (/kmCkm + h.C, 



2D 



i-D- j2 (/i/- 



T,(/m) 



+ -^ek(A, AA) sin [2^/:D^k(A, AA)] | i? (/L/k™)^^^ + {fUm) 



where the second terms of the r.h.s. of Eqs. (|ij29|) and liJ30|l have been neglected 
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APPENDIX C: TRANSFORMATION OF THE ONE-PARTICLE OPERATORS 

To determine the transformation of the one-particle operators we have again to apply the unitary transformation 
(|14|l which was used before to renormalize the Hamiltonian. As in Ref. 8 we first make the simplest operator ansatz 
for the A dependent c creation operator 

cL(A) = «kacL + «k,A/L (Cl) 

with the initial parameter values corresponding to the unrenormalized operators (A = A) 

Mk,(A=A) = 1, '"k,(A=A) = 0. (C2) 

Because the A dependent operators have to fulfill the same anticommutator relations as the unrenormalized operators, 
one concludes that 



1 - \uu.xf + D\vu,x\^ 



(C3) 



holds for all k and A values. Thus, the transformation of the / electron creation operator is given by 

/L(A) = -Dvu,xcl^ + u^,JL. (C4) 

Thereby, the approximation [/km'/km]+ ~ ^ "^^^ used. 

To derive renormalization equations for the parameters Wk,A and Wk,A of the one-particle operators we again consider 
the transformation step from A to (A — AA). As in the case of the Hamiltonian [compare Eqs. H31|) and H35|l ] we obtain 
two equations for cjj.„j(A — AA) 



cL(A-^-A) = ^ik,(A-AA)cL + «k,(A-AA)/L 



,/, , pXx^&x t „— ^A,AA I „,, , p^A.AA f\ „ — ^A,AA 

"k,A t. Cj^^ e -t- fk.A e y^m ^ 



^Xx ^x ?t c — Xx 



(C5) 

(C6) 

where the first one is derived from ansatz (|C4|) . The second equation follows from the application of the unitary 
transformation ifTB)) to cj^„(A). To calculate the transformed operators in Eq. (jC6|) one has to retrace the procedure 
of appendix IB] so that we obtain 



„-f A.AA „t „ — Xa.AA — „i 



L + ek(A, AA) {cos [VdA^{X, AA)] - l} 4 
^ek(A, AA) sin fV:DAk(A, AA)] /^„, 

\/L> L J 



pXxAX ft p — Xx.AX 



fl„^ + ek(A, AA) {cos [^A^{X, AA)] - l} fl 
- VB&^iX, AA) sin \VdA^{X, AA)] 4„. 



(C7) 



(C8) 
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Finally, inserting Eqs. (jC7|) and IjCSp into (|C6|) we find the renormalization equations for the parameters Mk,A and 
vux, 



■Wk,(A-AA) - Wk,A = 

= Uk.AJcos y/DAk{X,AX) - l| - Vl3t;k,Asin Vi^Ak(A,AA) 

Wk,(A-AA) - Wk.A = 

= Wk.AJcos /DAk(A,AA) -l| + ^ukasin \/i?Ak(A,AA) 



(C9) 
(CIO) 
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